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Abstract 

In this paper we explore the idea that the chiral phase transition in QCD can be 
described as a transition from a disordered instanton liquid to a strongly correlated phase 
of polarized instanton anti-instanton molecules. We calculate the degree of polarization of 
the molecules as a function of the temperature and show that the resulting T-dependence 
of the fermion determinant drives the chiral phase transition. We also show how the 
polarization of the molecules can lead to a non trivial behavior of the energy density and 
pressure. Finally, we study the effect of the presence of molecules on the propagation of 
quarks at at T ~ Tc. We derive the corresponding effective interaction and find that the 
strength in the scalar-pseudoscalar channel is four times the strength in the vector-axial- 
vector channel which agrees with recent lattice QCD simulations. We give results for the 
quark condensates as well as mesonic and baryonic correlation functions and find that 
the 'screening masses' of chiral partners become equal for T > Tc, where we still observe 
substantial attraction in the scalar-pseudoscalar meson channels. 
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1 Introduction 



Since the first suggestion (see [0] and references tlierein) tliat instantons are related to 
tlie breakdown of chiral symmetry in the QCD vacuum, significant efforts were made 
P, ^, ^, 1^ to transform this idea into a quantitative argument. 

Recently, two results have led to significant progress in connection with these efforts. 
First, it was shown 0| that even the simplest possible instanton-based model of the 
QCD vacuum, the so called Random Instanton Liquid Model (RILM), predicts correlation 
functions which agree both with phenomeno logical information (see the recent review ^j) 
and lattice calculations 0. In particular, it was shown that many hadrons (including, 
e.g. pions and nucleons) are bound by the instanton-induced interactions, and that their 
properties (masses, coupling constants and wave functions) are reproduced by the model. 

Secondly, important new results were obtained from the study of 'cooled' lattice con- 
figurations. 'Cooling' is a procedure that relaxes any given gauge field configuration to 
the closest 'classical component' of the QCD vacuum. As emphasized in earlier works, the 
resulting configurations were found to be of multi-instanton type [jlO|. The recent work 



by Chu et al. [[1^ now concludes that the typical instanton separation is given by i? ~ 1.1 
fm while the typical size is about p ^ 0.35 fm. These numbers essentially reproduce the 
key parameters of the 'instanton liquid' picture of the QCD vacuum, originally suggested 
on purely phenomenological grounds 0. 

In addition to that, it was found that correlation functions as well as hadronic wave 
functions in most channels remain practically unchanged after 'cooling'. This confirms 
that the agreement of previous lattice calculations with the instanton model was not ac- 
cidental. In fact, these works provide a decisive experiment as far as the validity of the 
instanton model is concerned, demonstrating that even after eliminating such QCD phe- 
nomena as the perturbative gluon exchange and confinement, one still observes hadronic 
bound states. Moreover, the corresponding masses and wave functions appear to be only 
mildly affected! 

Although many details remain to be worked out, it seems fair to say that instantons 



are indeed a major component of the QCD vacuum, producing the quark condensate, the 
low lying hadronic states and many other nonperturbative features of QCD. 

In this paper we discuss properties of the instanton ensemble at finite temperature. 
Our focus is especially on the region around the chiral phase transition T ~ T^. The first 
attempt to understand this phase transition as a rearrangement of the instanton liquid, 
going from a random phase at low temperatures to a strongly correlated 'molecular' phase 
at high temperatures was made in |T2|. The essential point is that while individual instan- 
tons are strongly suppressed in the chirally restored phase (because of the corresponding 
zero modes), instanton— anti- instanton (//) pairs have a finite probability even above Tc. 
This idea was recently reexamined in [|T^, and although the basic ingredients of the model 
are the same, the two works suggest significantly different scenarios for the chiral phase 
transition. 

According to the original idea, chiral symmetry is restored because most of the in- 



stantons are removed by a finite temperature instanton suppression factor ||I4|, |15|. This 
factor describes the effect of Debye screening on the fluctuations around the instanton 
solution. Of course, it also affects the high temperature molecular phase and allows only 
a small number of molecules above T^. In the more recent work |T^, on the contrary, 
it was argued that screening should not be a dominant effect at temperatures below or 
just above Tc (see also |TB[)- Instead, it was shown that the temperature dependence of 
the fermion determinant can provide chiral restoration even without any additional sup- 
pression factors^. This means that one may have a significant number of instantons even 
above Tc, and that the cause of the phase transition is not the suppression of instantons 
but a rearrangement of the instanton liquid. 

This new scenario for the chiral phase transition implies a number of nonperturbative 
effects in the region just above Tc. One of these effects was already mentioned in 



instantons can provide a substantial contribution to the energy density and pressure 
even above the chiral phase transition. In this paper, we want to study the effects of 



^In earlier works on the subject [|2[ |l^, |l^ the high temperature instanton suppression factor was 
extrapolated to all temperatures. 



the rearrangement of the instanton vacuum in more detail. In particular, we consider 
the behavior of condensates and hadronic correlation functions through the chiral phase 
transition. For this purpose we consider a schematic model in which the instanton liquid 
is described as a mixture of a random and a molecular component. The phase transition 
is then studied as a function of the fraction of correlated instantons. 

The paper is organized as follows. In section 2 we give a general discussion on cor- 
relations in the instanton liquid and introduce the concept of instanton— anti-instanton 
molecules. In section 3 we study the phenomenon of polarization of molecules at finite 
temperature and in section 4 we describe how this effect determines the energy density 
and pressure of the instanton ensemble. A general treatment of the new interaction in- 
duced by the presence of molecules is given in section 6 while in section 7 we present the 
results of a simulation of a mixed ensemble of random instantons and molecules. 

2 Correlations in the instanton ensemble 

Let us introduce the necessary general formulae and notations. The ensemble of interact- 
ing instantons is described by a partition function of the type 

^ = E WJWl / n id^^ diP^)pf'] eM-S^nt) U det(D + mf) (1) 

N+N- -'V + '-'V-- •' i f 

where A^_|_ and are the numbers of instantons and anti-instantons, dfli is the measure 
in the space of collective coordinates (12 per instanton in SU{3)) and 

4p) = C„„p-(iJ)^"%xp(-i^) (2) 
_ 4.6exp(-1.867V,) 

is the semiclassical amplitude for a single instanton which depends on the running coupling 
constant g{p)- Furthermore, Sint is the gluon-induced interaction, which we do not specify 
in this work. Instead, we focus on the last factor, which appears after the integration over 
fermions has been carried out. As is the case for the bosonic zero modes, the integral 
over the fermion zero modes has to be performed exactly. Assuming that the fermion 



determinant factorizes in a zero mode and a nonzero mode factor, the first factor is given 
by det(TT"l' + mj), whereas the second factor is taken into account to gaussian order. The 
matrix Tjj is the N/2 x N/2 'hopping' matrix 

Tjj = j d^x (l)\^{x - zj)ib^(t)io{x - zi). (4) 

Here, 2;/, zj are the centers of the instantons and for simphcity we assume thatg, = 
= N/2. The statistical mechanics described by the partition function (|I]) is comph- 
cated and so far direct simulations have only been carried out at zero temperature. 

The instanton solutions and the corresponding zero modes are known analytically for 



non-zero temperature |T9[, and detailed studies of the temperature dependence of the 
'hopping' matrix elements were performed in The general structure of these 

matrix elements is given by Tjj = U4/1 + {uf/r)f2. Here, the fourvector = (u, U4) is 
defined by Ujj = m^t^, where Ujj is the upper 2x2 corner of the N^ x N^. matrix that 
describes the relative orientation of the instanton and the anti- instanton and = (r, —i). 
The two invariant functions /i and /2 can be parameterized in the form 

. ttT sin (vrrT) cosh (vrrT) 

/i = «7 1 / ^ ■Fi{T,r,T), (5) 

(cosh (vrrT) — cos [tttT) + (T))^ 

. ttT cos (vrrT) sinh (vrrT) 

(cosh (vrrT) — cos {tctT) + K2(T))^ 

where r and r are the separation between the instanton and the anti-instanton in the 
spatial and time direction. All quantities are given in units of the geometric mean of 
the two instanton radii ^JpipA- The functions Fi, F2 and ki, K2 provide the correct 
normalization in the limits of zero and infinite temperature. They are explicitly specified 
20[| and used in our simulations but their detailed structure is not important for the 



m 



purpose of our arguments. 

Neglecting all interactions one arrives at a random ensemble of instantons in which the 

distribution over all collective coordinates such as positions and orientation is given by 

the corresponding invariant measure. As was mentioned in the introduction, this model 

leads to a simple and phenomenologically successful description of the QCD vacuum. 
^Recall that any condition applied to a macroscopically large system should not be important anyway. 



Notable exceptions are those channels in which instantons produce a strong repulsion, in 
particular the rl and 5 (scalar-isovector) channels. 

Correlations in the instanton liquid have two aspects, which we will refer to as long 
and short range correlations. Long range correlations are related to the phenomenon of 
screening of the topological charge. Its most important manifestation is the vanishing of 
the global topological charge (and the vanishing of the topological susceptibility), which 
takes place provided the theory has massless fermions and one is considering a sufficiently 
large volume. As expressed by the Witten-Veneziano relation, these correlations are re- 
lated to the physics of 77' meson. We plan to discuss this question in a separate publication 



In the present work we want to focus on short range correlations, which can be stud- 
ied by considering 11 pairs. Correlations in the instanton liquid that may lead to the 
formation of II pairs are caused by the gluonic action Sint or the effective fermion action 



flavors A'^^ in the theory. However, both types of interaction favor the same relative orien- 
tation of the instanton and the anti-instanton. This orientation is most simply described 
in terms of an angle Q defined by 



where = (r, r) is the vector connecting the centers of the instantons and is the 
relative orientation vector introduced above. Maximal attraction, for both the gluon and 
the fermion interaction, corresponds to ^ = 0. 

Since the probability of a given configuration is proportional to the value of the fermion 
determinant, configurations with a large determinant should be preferred. For a single 
II pair we have det(D) ~ {cosO)'^^^ . If this // pair sits in the vicinity of = we will 
refer to it as a 'molecule'. Clearly, it maximizes the fermion determinant, or minimizes 
the corresponding 'energy', — log(detZ)). 

However, dealing with a statistical system one should not simply minimize the energy, 
but rather take into account the competition between minimizing the energy and maxi- 



21- 



Sf = A^/Tr log(Z)). Naturally, their relative role depends strongly on the number of light 



cos(^) 



u ■ R 



(7) 
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mizing the entropy (i.e. minimize the free energy). In other words, we have to determine 
the relative weight of the most attractive (molecular) orientation relative to all others in 
the SU (3) group from summing over all states in the partition sum. 

Before we consider this problem, let us make a brief digression. Even if the contribu- 
tions of 'molecules' in the QCD vacuum is not large, there are external parameters that 
can increase their role. One way of doing this is to consider QCD-like theories with an 
increasing number of light flavors Nf. As the fermion determinant is raised to a higher 
power Nf, the role of correlations induced by the determinant certainly increases. Thus, 
one may anticipate the existence of some critical number of flavors Nf^ cr above which the 
instanton ensemble is dominated by molecules ||2^, In this paper, however, we want 



to consider another external parameter that has a significant effect on correlations in the 
instanton liquid, namely the temperature. 

3 Polarization of the // molecules 



In reference |]T3[ a schematic model was developed in order to deal with the complicated 
statistical mechanics described by the partition function (|Tp. In this model the instan- 
ton ensemble is described as a mixture of a molecular and a random component. The 
differential activities for the two components are assumed to be 



dZ^ = C''dp,dp2d''RdUip^p2f-'exp[-K{pi + pl){plna + 2plnm)\ ((T,jT;,)^0 (8) 

for the molecular component and 

dZ, = 2Cdpd''RdUp'-'exp[-Kp\'^,na + 2^nj] {j^TiTTYf/' (9) 

for the atomic component. Here, n^, denote the densities of the random and the molec- 
ular components, p^, are the average square radii of instantons in the two components, 
C is the normalization of the single instanton density (see (^) and b = y-^c ~ |A^/ is the 
coefficient of the Gell-Mann-Low function. 

The model uses a simplified gluonic interaction corresponding to an average repul- 
sion (Sint) = i^PiPi parameterized in terms of a single dimensionless constant n. The 



fermion determinant for the random component is approximated by (-^rTrTT^)^-'/^ where 
the average is over all positions and orientations. For the molecular component, on the 
other hand, the overlap matrix element is first raised to the Nf power and then aver- 
aged {{TjjTjj^ ^). This average is only performed over positions, whereas the relative 
orientation is kept fixed. 

Using the overlap matrix elements defined in the last section, one finds a very specific 
temperature dependence of the quark induced interaction. The average determinant for 
the random component gradually decreases with temperature, whereas the determinant 
for the molecular component first increases (at T ~ T^) and eventually, at larger T, 
starts to decrease. This behavior leads to the disappearance of random instantons at 
high temperatures and to a phase transition to a purely molecular phase in which chiral 
symmetry restored. 

As we are going to show in this section, the reason for this temperature dependence is 
given by a strong and rapid polarization of the molecules in the critical region. Polarization 
in this context means that at finite T the vector connecting the centers of the instanton 
and anti-instanton = zj^ — zj^ = (r, r) is no longer distributed isotropically (as it is 
at T = 0), but tends to point in the time direction. Roughly speaking one may say that 
at T ~ Tc the "instanton liquid" actually becomes a "liquid crystal" of the nematic type. 

In general, the anisotropy of a field theory at finite temperature is of course a conse- 
quence of the different boundary conditions imposed in the time and spatial directions. 
An important consequence of these boundary conditions is the fact that the propagators 
for bosons and fermions become qualitatively different at finite temperatures. Both are 
anisotropic, but fermions are subject to a qualitatively new phenomenon, the exponen- 
tial screening of the propagator in the spatial direction. For a free massless fermion, the 
propagator is S{r) ~ exp(— vrTr) where vrT is the lowest Matsubara frequency. There- 
fore, with increasing temperatures quarks develop a strong preference to propagate in the 
time direction. This is the main mechanism that produces an anisotropy of the instanton 

ensemble at finite temperature 0. 

•^The gluonic //interaction at finite T was studied in ||20|] , with the result that it remains approximately 
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Before we come to quantitative results, let us give a qualitative explanation why the 
phenomenon takes place. Suppose the vector connecting the I and I centers is = (r , r) . 
Using the overlap matrix elements specified above, the probability to have such a molecule 
is roughly proportional to 

det(iD) ~ |sin(7rrr)/cosh(7rTr)p^^. (10) 

This factor is maximal^ for r = and t = (3/2 = (1/2T), which is the most symmetric 
position of the II pair on the torus. Note that it corresponds to a molecule polarized in 
the time direction. As a function of temperature one may expect the strongest effects to 
occur when the temperature is such that the period of the torus in the direction is about 
twice the typical size of a molecule. In the euclidean formalism, the phase transition 
becomes a simple geometric effect: the transition occurs when the // molecules optimally 
fit into the 'temperature box'. 

In order to study this effect in more detail we have performed a simulation of the 
partition function (|l]) for a single II pair at finite temperature T. The interaction is given 
by the fermion determinant discussed in section 2 and the gluonic interaction specified in 
|[20|1 . The simulation is based on a simple Metropolis algorithm in which a new random 
configuration of the II pair is accepted with a probability p ~ exp{—Seff), where the 
effective action is the sum of the gluonic and fermionic parts, S^jf = Smt + log(detZ)). 
Since there is only one pair present in our simulation we keep the size of the instantons 
fixed while the positions and orientations are updated. 

In order to characterize the orientation of the II molecule we introduce the angle a 
between the vector (connecting the / and / centers) and the (3-d) spacelike plane: 
tana = i?o/|-R| = t/t. The distribution of the polarization is then determined by av- 
eraging the a over a large number of configurations. In fig. 1 we show the differential 
distribution of the angle a for different temperatures and two different cases: (i) two light 

and one heavy flavor (QCD); and (ii) four light flavors. The distributions are normalized 

isotropic until rather high temperatures. 

^Note also, that near the maximum det(il)) ^ exp[— iV/7r^r^((5r^ +(5r^)], so at large number of flavors 
one can use a saddle point method to evaluate the contribution of this region. 
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to the isotropic 4-d distribution, which means they are divided by the relevant Jacobian 
cos^ a. In both cases we observe that the distribution is practically flat (=isotropic) at 
low temperatures T < 100 MeV, but becomes strongly peaked at a = 7t/2 for larger T. 
At very large temperatures, this peak starts to disappear: this happens because at high 
T the time dimension is so strongly squeezed that no molecules can be oriented this way. 
As expected, the polarization is stronger in the case of four light flavors. 

In fig. 2 we show fraction of molecules in the ensemble with \R\ < min(i?0; 1/T — Rq). 
We observe that the degree of polarization jumps to a very large value ~ 70% over a fairly 
small temperature interval T ~ (120 — 150) MeV. 

Finally, let us make a comment on the physical meaning of the polarization phe- 
nomenon in a less technical language. The // molecules are virtual (or failed) tunneling 
events, in which the gauge fields penetrate into a new classical vacuum only for a short 
period of time, and then return back. For that reason, they do not contribute to the quark 
condensate and other related quantities. "Polarization" of the II molecules at T ~ 
means that at such temperatures the tunneling is concentrated in the vicinity of the same 
spatial point. 

4 Magnetic versus electric fields, and the equation of 
state 

Lattice data on the thermodynamics of the chiral phase transition suggest that the pres- 
sure p{T) remains relatively small around the transition region, while the energy den- 
sity e(T) changes by about an order of magnitude over a small| interval 6T <^ Tc. A 



parametrization of the data obtained by Kogut et al. [|2^ together with a general dis- 
cussion can be found in |2^. These authors stress that at chiral restoration T ~ Tc not 
only the (relatively small) quark related energy and pressure are significantly changed, 
but that the (much larger) gluonic energy and pressure should also be modified. 



^ For three and more massless quarks the transition is beUeved to be of first order, but the phenomenon 
mentioned above is also observed for two massless quarks (when there is a second order phase transition) 
and for light but massive quarks, when there is no strict phase transition at all. 



Let us mention one more point emphasized in |2^: there seems to be a contradiction 
between the small transition temperatures ~ 140 MeV seen in unquenched lattice 
simulations and the large value of the phenomenological bag pressure B ~ 500MeV/fm^ 
derived from the trace anomaly 0, ^ (at T=0). With the usual bag model equation 
of state, in which the pressure at T > is written as the sum of the Stefan Boltzmann 
contribution for an ideal quark-gluon gas minus the nonperturbative bag pressure, one 
is unable to get a positive pressure at such at low transition temperature! The natural 
alternative is to assume that the nonperturbative phenomena responsible for the bag 
pressure at T = are partially present at T > as well, so that only part of this large 
B is removed across the transition region. 

The question whether the rearrangement of the instanton ensemble into a molecular 



phase can explain this behavior was studied in , where the schematic partition function 
defined by was used to determine the instanton related contribution to the equation 
of state. It was indeed found that instantons give a sizeable contribution to the energy 
density and pressure, that the two quantities behave quite differently and that the results 
are roughly consistent with the lattice data. In particular, it was shown that, (i) Pinst{T) 
jumps down, to about half of its value at zero temperature, while, (ii) emstiT) jumps up, 
thus contributing to the "latent heat" of the transition. 

The reason for this behavior was attributed to the specific temperature dependence 
of the fermion interaction used in that model. In the last section we have explained that 
this T-dependence is caused by the polarization of // molecules. Now we want to show, 
that this phenomenon leads to a simple microscopic explanation of the energy density 
and pressure, related to the (color) electric and magnetic fields associated with polarized 
molecules. 

Let us first recall some well known facts about instantons at T = 0. The instanton 
(anti-instanton) fields are selfdual (or anti-selfdual), Gfj_y = ±je^uaf3Gai3- In Minkowski 
space they correspond to classically forbidden paths and the electric and magnetic fields 
are related by Ek{x) = ±iBk{x). As a consequence, in the classical energy density ed = 

10 



{1/2){E'^ + B"^) the negative electric term is exactly compensated by the positive magnetic 
one, so the configuration may be created 'out of nothing', in agreement with the classical 
equations of motion. 

Furthermore, instantons contribute to the energy density of the QCD vacuum, but 
this happens not at the classical but at the quantum (one-loop) level. This is quite 
natural, since instantons correspond to tunneling phenomena, which are known to lower 
the ground state energy. In QCD, this fact can be expressed in terms of the 'stress tensor 
anomaly' 

which is non-zero for the instanton solution. Here, b = {ll/3)Nc — (2/3) A^/ is the first 
coefficient of the beta function. For a dilute ensemble, one can directly relate the energy 
density to the instanton density e = —\ninst [0- As we already emphasized above, this 
energy is numerically rather large. 

At finite instanton density (and especially in strongly correlated molecules) the fields 
are not purely self-dual. However, at T = molecules are unpolarized: therefore only 
Lorentz scalars can have non zero vacuum expectation values. Since the liquid is relatively 
dilute, these deviations from self-duality just lead to a relatively small correction to the 
right hand side of the anomaly equation. 

At finite temperature, instantons contribute to the equation of state in two ways. 
First, the instanton density will change, leading to a change in the bag pressure. Second, 
the polarization of // molecules in the time direction causes the expectation values of the 
electric and magnetic fields to be different. While the first effect gives equal but oppo- 
site contributions to the energy density and pressure, non-selfdual fields give a positive 
contribution to both the energy density and the pressure. 

In order to study this effect, we have calculated the color fields for a single molecule 
with given polarization. Since // configurations are not exact solutions of the euclidean 
equations of motion the quantitative result depends on the specific ansatz for the gauge 
fields. Here we show the two most widely used forms of the gauge potential, the sum 
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ansatz (full lines in figs. 3 and 4) and the ratio ansatz (dashed lines in figs. 3 and 
4). 

In fig. 3 we show the modification of the action Sint = Stot — 25*0 and the energy for 
a single molecule polarized in the time direction at a temperature of 200 MeV. Results 
are given for the most attractive and the most repulsive orientation. Indeed, one finds 
considerable attraction for the favored orientation cos 6 = 1. This interaction is weaker 
in the ratio ansatz, but the qualitative features are the same as in the sum ansatz. The 
new result is the lower panel, showing the polarization of the fields in the euclidean 
time direction. Both in the sum and in the ratio ansatz cases, the fields are dominantly 
magnetic and give a positive contribution to the classical euclidean energy density e = 
l/2((i?^ — E"^)). The effect is on the order of 15% in the sum ansatz, whereas in the ratio 
ansatz it is roughly 10%. 

In fig. 4 we show the same calculation for a molecule which is oriented along a spacelike 
direction. Again we observe attraction for cos^ = 1, but this attraction is considerably 
weaker as compared to a polarized molecule. Also, in the ratio ansatz the molecule is now 
essentially self dual, while in the sum ansatz it is actually dominantly electric. 

Let us now come back to the discussion of the energy density and pressure 

e = 1(B^-E')-9'y^,{E' + B'), (12) 

P = \{B^-E^) + g^-^^{E^ + B'). (13) 
1287r"' 

The important point is that although {E'^ + B'^) is much larger than {B"^ — E"^), the former 

quantity is suppressed by a smaller coefficient. Therefore, even a moderate deviation from 

selfduality leads to a substantial contribution to the energy density and pressure! Taking 

{B"^ — E'^)/{B'^ + E"^) ~ 10% from the calculation of a single molecule, one can get 

e ^ 225MeV/fm^ even if the density of instantons does not change across the transition. 

Thus it appears that the observed jump in the energy density can be explained by a change 

in the (instanton related) bag pressure plus the contribution from instanton interactions 

and the perturbative Stefan Boltzmann part. 
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5 Quark interactions induced by // molecules 



Instantons give a very successful phenomenology of the QCD vacuum, because they pro- 
vide a mechanism for chiral symmetry breaking and generate strong (and quite specific) 
interactions between light quarks. Although instantons do not lead to confinement the 
interactions are able to bind quarks into mesons and baryons and lead to a quantitative 
description of the observed spin splittings. At the one-instanton level the interaction is 



described by the famous 't Hooft effective interaction |28 



^ ~ l[i^<l>o){M) (14) 
/ 

where are quark operators and 0o is the zero mode wave function. The product runs 
over the light quark flavors. This is still a complicated nonlocal interaction which depends 
on the color orientation of the instanton. After taking the long wavelength limit and 
averaging over the color orientation it can be reduced to a local 2A^j-point interaction 
3^, 3D|. In the case of two light flavors it reads 



where n{p) denotes the density of instantons. Here, ip is an isodoublet of quark flelds 
and the four vector r~ has components (r, i) with r equal to the Pauli matrices acting in 
isospace. The Dirac matrix a^u is deflned by (t^,^ = ^[y^, ju]- The interaction ([T5|) success- 
fully explains many properties of the (T=0) QCD correlation functions, most importantly 
the strong attraction seen in the pion channel. 

At finite temperature the collective coordinates of instantons are no longer random, 
but become correlated. Therefore, it is no longer possible to average over the instanton 
orientations and use the effective Lagrangian (^5|). However, due to the presence of 
molecule there should be strong non-perturbative effects even around Tc s. For this 
reason we would like to investigate what kind of interactions are induced by correlated 
// pairs. We first study the quark propagator in the field of an instanton anti-instanton 
molecule and then derive the effective vertex that reproduces all four quark correlations. 

13 



The quark propagator in the field of a single molecule is given by 

5j (x, y) = So{x, y) + ^_ {riM^'^v) + • (16) 

Here 0/ and denote the zero mode wavefunctions in the field of the instanton and the 
anti-instanton and Tjj is the corresponding overlap matrix element. For simplicity we use 
the zero temperature profiles. We neglect all non-zeromode contributions except for the 
free propagator Sq. The propagator (p!6D satisfies the following symmetry relation 

S\x,y) = S{y,x). (17) 

For Nc = 2 the symmetry group is larger. In that case we can relate the propagator to 
its transpose by 

Six,y) = -a2'y5CS'^{y,x)C-f5a2 (18) 

which is also known as the Pauli-Giirsey symmetry [^. The Pauli matrix a2 acts in 
color space. A remnant of this symmetry is still present for more than two colors. If 
we remember that correlators of color singlet currents only depend on the relative // 
orientation, which in the case of molecules is given by R^t^ , it is clear that the symmetry 
(|l^) holds in the gauge that one of the pseudoparticles has the identity as orientation 
matrix. In that case ct2 contains a projector of SU (Nc) on the SU (2) subgroup where the 
// molecule has its support. 

To investigate the consequences of the symmetries 
decomposition of the propagator 

S^\x,y) = Y.afix,y)T„ (19) 

i 

where Fj is the complete set of (hermitean) Dirac matrices. Due to the chiral structure 
of the propagator only the vector ay ^ and the axial vector a^^ components are non-zero. 
For the scalar, pseudoscalar and tensor components we have 

as = ap = aT^^u = 0. (20) 



17, 18) let us consider the tensor 
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As a consequence the scalar and the pseudo-scalar correlator are degenerate. Also 
the axial and the vector correlator are degenerate, whereas the tensor correlator is not 
affected by the presence of molecules. 

The symmetry (|17D leads to the relation 

ai{x,y) = al{y,x). (21) 

Under the restrictions quoted below (^) also the following two equalities hold 

avt,{x,y) = -(T2a^^{y,x)a2, (22) 
aAi,{x,y) = a2a^^{y,x)a2, (23) 

An immediate consequence of ( p^ ) is that 

Tt{S{x,x)j^) = 0. (24) 

This term enters the disconnected contribution to the isoscalar correlation function. The 
relation (l2^) therefore implies that in the molecular vacuum the rho and omega mesons are 
degenerate^]. To exploit the relations (^,^) further, we remind the reader the definition 
of a flavor non-singlet meson correlator 

n*^(a;, y) = Tt{S{x, y)TS{y, x)T), (25) 

and the corresponding diquark correlator 

n?(x, y) = Tt{S{x, y)TCa2S''{x, y)a2CT), (26) 

where we have chosen a particular color component of the diquark (every antisymmetric 
color matrix defines a diquark current). Again, a2 projects onto a SU{2) subgroup of 
SU{Nc). From (|T^) we conclude that meson correlators in the channel F are degenerate 
with diquark correlators in the channel F75 up to a color factor. This factor is given 

by 2/Nc{Nc — 1) which counts the total number of embedings of instantons in SU{Nc). 

^This is quite remarkable since the near-degeneracy of the rho and omega mesons has no really good 
explanation. In fact, finite-T QCD sum rules suggest that the degeneracy is lifted at finite T whereas 
our model predicts that it remains present. 
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For Nc = 2 the Pauli-Giirsey symmetry becomes exact and meson and diquarks are 
degenerate. 

In it was shown that all meson correlators can be expressed in the coefficients of 
the tensor decomposition of the propagator. In the molecular vacuum, the coefficients 
obey simple relations which allow us to find relations between the correlators. We find 

El4V(^,i/)P = 3EK(^,i/)^ (27) 

kl kl 

and 

kl ^ kl 

where the averaging is performed over the color and spatial orientation of the molecule 
and we have assumed that x and y are far away from the locations of the pseudoparticles. 
The relation (^71) determines unflavored correlation functions P], whereas (p8|) gives the 
disconnected contribution for flavor singlet correlation functions. Note that (|2^) implies 
that ay^{x,x) = 0. If the propagator is calculated in the spatial direction, the relations 
(^,^) remain valid even the molecules are polarized in the temporal direction. 

The effective vertex for quark-quark scattering in the molecular vacuum can be ob- 
tained using the standard reduction formula 

<PlP2\Cmol\p3PA> = - I nrfV,e^Pl^'l+*^'^^'^-*^^^^-^^'^^^C(j3i)™V^^(p2)n, 

< 0|T(g7(xi)g5(x2)g^(a;3)gf(x4))|0 > {p,)kp%{PA)i,^i, (29) 

where ip^ is an isodoublet fermion field with color index a and spinor index i. If we treat 
the correlated instanton ensemble as a dilute gas of molecules, the vacuum expectation 
value can be calculated using the quark propagator (|16|). We only have to perform the 
average over the global color orientation of a molecule whereas the relative orientation 
between the the instanton and the anti-instanton is kept fixed. At low temperatures, 
molecules are unpolarized and we average over the direction of the vector connecting the 
centers of the instantons. After taking the long wavelength limit we obtain a local four 
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fermion interaction 



mol 



G 



Nc{N, - 1) 

"iV,(iV2-l) 
1 



iVc-2 
Nc{N^ - 1) 



iVc(iVe - 1) 

with the couphng constant 
G = 



n{pup2) dpidp2 -i2-(27rpi)2(27rp2) 



ii^ii^i^Y - (^^7^75^^)^ 



(30) 



(31) 



Here n{pi,p2) is the total tunnehng probabihty for the // pair. Note that the isospin 
structure of ([3^ ) is different from ([T5|), is a fourvector with components (r, 1). When 
calculating correlation functions from this effective vertex both the direct and the ex- 
change term have to be taken into account. In order to extract the relevant coupling in 
each channel directly from the Lagrangian we add the exchange term and Fierz rearrange 
it into an effective direct interaction. The resulting Fierz symmetric Lagrangian reads 

2 



mol sym 



1 



2m 



{^T'^^^^Pf + {^T^l^^l^i^f] + ^(^^7^75^^)' + A 



(32) 



G 



where £3 denotes the color octet part of the interaction 

iV,-2 
2iV,(iV2-l) 
1 

4iV,(iV, - 1) 
1 

2iV,(iVe - 1) 

2Nr - 1 



+ 



(^r»A*7M^)' + (^r'^A*7M75^)' 
mi.i^f + (^A^7m75V')' 
{i^^l.i^f - (V^A^7m75^)'" 



(33) 



2iV,(iV2 - 1) 

At temperatures around the critical temperature for the chiral phase transition, molecules 
become polarized in the time direction and we should not average over the relative 
position of the instanton and the anti-instanton. In this case we get an effective la- 



grangian which is identical to (30,33) but with all vector interactions modified according 



to (^/'7m^)^ ^ 4(^-70^)^ 
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Like the effective interaction in the dilute instanton gas, (p!5[) , the result for the molec- 
ular vacuum is SU (2) x SU (2) symmetric. But in addition to being chirally symmetric, 
the effective lagrangian in the molecular vacuum is also U{1)a symmetric. Local fermion 
lagrangians of this type have been studied extensively in the context of the Nambu and 
Jona-Lasinio (NJL) model [Q. However, in contrast to the NJL model the effective in- 



teraction in the instanton model changes above the chiral phase transition. Furthermore, 
the NJL model has four unknown coupling constants, even when one restricts oneself to 
zero temperature and colour singlet terms only. At finite temperature, the number of 
independent couplings is even larger since the space- and timelike vector currents can 
have different couplings constants. In our model, all coupling constants are specified in 
terms of a single parameter G. 

Let us now compare the 't Hooft interaction ( [T3| ) (dominant at T = 0) and the 
'molecule-induced' one we just derived, valid for T > T^. The main effects of the 't Hooft 
interaction are (i) the strong attraction in the pion channel and (ii) the strong repulsion 
in the t]' channel. The 'molecule-induced' effective lagrangian is also attractive in the 
pion channel, but (since chiral and U{1)a symmetry are restored) the same attraction 
also operates in the scalar-isoscalar (a), scalar-isovector (5) and pseudoscalar-isoscalar 
[r]') channels. 

Furthermore, this lagrangian also includes an attractive interaction in the vector and 
axial vector channels. If molecules are unpolarized, the corresponding coupling constant 
is a factor of four smaller than the scalar coupling. If they are completely polarized, the 
attraction is present only in the operators containing 70, but becomes equally strong as 
the attraction in the scalar-pseudoscalar channels. Transversely polarized vectors are not 
affected in this case. 

There are further non-trivial qualitative features of the new effective Lagrangian. Note 
that it predicts no splitting between the isoscalar (u) and isovector (p) vector channels. As 
we explained above, this is a consequence of the remnant of the Pauli-Giirsey symmetry 
in the molecular vacuum. At the same time, there appears additional repulsion in the 



axialvector-isoscalar (/i) channel compared to axialvector-isovector (ai) case. 

In order to estimate the effective interaction in baryons we would like to start with 
a simpler problem and study the interaction for diquarks. The effective lagrangian (|30|) 
can be rearranged into the form 

2 



L = G 



- {i^l^l^Cr^ii'^r ) C^,^^T^(i'^i,)\ I + £6 (34) 

where L% denotes the effective interaction among color sextet diquarks. Only the anti- 
symmetric color matrices enter in (^4|). In order to facilitate the comparison with the 
mesonic interaction ||3^, they are normalized according to Tr(/3^/3^') = N^.S'^^' . Using this 



normalization, the strength of the interaction as compared to the free correlation function 
can be inferred directly by comparing the coupling constants. For diquarks, the isospin 
wavefunction is determined by the dirac structure @, |33| . Here, denote the symmetric 
Pauli matrices and = is the antisymmetric one. 

From the effective lagrangian ( ^If ) we find an attractive interaction for scalar ip'^C'y^ip 
and pseudoscalar ip'^Cip diquarks. The interaction for vector diquarks is also attractive 
with a coupling constant that is a factor 4 smaller than the one in the scalar channel. All 
coupling constants are a factor (Nc — l) smaller than the ones in the corresponding mesonic 
channel. Again, this is a consequence of the remnant of the Pauli-Giirsey symmetry. 

The nucleon can be thought of as consisting of an equal mixture of scalar and vector 
diquarks coupled to a third quark, whereas the delta resonance only contains vector 
diquarks. Using the effective interaction derived above, we would therefore expect an 
attractive interaction for the nucleon. The delta-nucleon splitting, however, is expected 
to be significantly smaller than the pi-rho splitting since the coupling in the diquark 
channel is a factor of 2 smaller than the coupling for mesons. 

The most direct consequence of the Lagrangian (|30D is the behavior of mesonic cor- 
relation functions. We will study this question in the next section. The corresponding 
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measurements on the lattice provide the relevant test in order to check the correctness 
of our model. Before we come to this, we would like to consider an even simpler set of 
physical quantities at nonzero T, the expectation values of four quark operators. These 
matrix elements play an important role in QCD sum rule calculations and some of them 
can be considered as additional order parameters for chiral symmetry breaking. It is 
therefore interesting to check whether they vanish at T = in a universal manner. At 
low temperatures, the temperature dependence of these expectation values can be de- 
termined by calculating the contribution from single pion states. At large temperatures 
T ~ Tc, however, little is known about the T-dependence. Usually it is assumed that the 
expectation values can be factorized 

^^^^^^^'^^^ ^ (mlNfy ^^'^^^^'^^'^ ~ ^'^^^'^^ ^^^^ 

where the double brackets indicate thermal expectation values. This would imply that 
all expectation values of four quark operators are restored at T^, and that they vanish at 
twice the rate of the quark condensate. 

Alternative SU{Nf) order parameters are given by the expectation values 

Oi(T) = {{[uL^QUL-dL^odL][L^ R])), (36) 

02(T) = {{[uL-fiUL-dLl^dL][L^R])), (37) 

0,{T) = {{[uLiot^UL~dLlot''dL][L^R])), (38) 
0,{T) = {{[uLlit''uL-dLiiedL][L^R])), (39) 

where = A"/2 are the generators of SU{Nc). These operators enter the Weinberg sum 
rules at non-zero temperatures, related to the difference between vector and axial corre- 



lators Instead of the operators listed above, one may also consider the combinations 
only containing left or right handed quark fields. These operators are related to the sum 
of the vector and axial vector correlators, and need not be restored above the chiral phase 
transition. 

A quantity sensitive to U{1)a violation is nothing else but the expectation value of 
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the 't Hooft vertex 



O'tHooft^j.^ = ((det(g7,^g;,,))), (40) 

where the determinant is taken over the different quark; ffavors f = u,d, . . . and L, R stand 
for left and right components of the quark fields. Around or above there is no quark 
condensate, and 'unpaired' instantons have very small density 0{m^f) [^]. However, in 
the measurement of quantities like this one, the denominators of the quark propagators 
cancel the current quark masses, recovering the famous 't Hooft effective interaction. In 
a sense, the operator considered can induce a tunneling event by itself. Unfortunately, 
this implies that a small fraction of the configurations will produce a large signal. This is 
certainly not an easy way for measurements to be performed. 

Finally, at the end of this section let us address the question of the absolute strength of 
the 'molecule-induced' new effective interactions. Instead of 7 unknown constants (which 
would appear in a general NJL-type Lagrangian for T > T^) we have only G, which is 
determined by the number of // molecules nmoi{T). 

Let us make a few general remarks. First, there should be a general upper bound, 
nmoi{T) < ric (the density should be smaller than some critical number for all T > T^) 
following from the condition that the effective interaction should not be strong enough 
to lead to a rearrangement of the vacuum and cause spontaneous breakdown of chiral 
symmetry. Second, in a fully polarized configuration the interaction in the pion and 
longitudinal vector channel are identical. This might lead to a situation in which all 
scalars as well as the longitudinal components of vector and axial mesons are massless. 



This speculative scenario resembles the "vector limit" studied in |]35[ . 

Third, and the most natural case, is that the density of molecules and the effective 
interaction is simply too weak to produce any bound state. However, it can still be 
important and lead to observable effects: one of them being modification of the 'screening 
masses' to be discussed in the next section. 
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6 Correlation functions in the 'Cocktail Model' 



Having outlined the main qualitative features of the interactions induced by // molecules, 
and local averages induced by them, we now want to study more complicated predictions 
of this model involving various correlation functions. 



For this purpose we introduce a schematic model along the lines proposed in |T3 
In this model we assume that the total number of instantons at temperatures around 
the chiral phase transition is not much smaller than at T = 0. We then study the 
phase transition by varying the fraction / = 2Nmoi/Ninst of instantons correlated in 1 1 
molecules. 

We have carried out our simulations in a box 5.8'^ x 1.3 fm^, corresponding to a total 
volume V = 256 fm^ and a temperature T = 150 MeV. This is the temperature that 
determines the boundary conditions on the fields and enters in the calculation of the 
propagators. As we have explained in section 4 we consider a scenario in which the phase 
transition takes place over a fairly narrow interval of temperatures 6T <C Tc. Since the 
instanton fields vary smoothly with temperature, we keep the size of the box fixed while 
we vary the fraction of molecules. 

The calculation of the spectrum of the Dirac operator, the quark condensates and 
hadronic correlation functions follows our earlier work on the random instanton model 
1^. The essential idea is to exactly diagonalize the Dirac operator in the space spanned 
by the zero modes of the individual instantons and treat the effects of the non zero modes 
on the single instanton level. In this work we simply replace all the eigenfunctions and 
overlap matrix elements by their finite temperature counterparts. This means in particular 
that the quark propagator is antiperiodic in the time direction and becomes screened for 
large spacelike separations. 

In fig. 5 we show the behavior of the quark condensate {{qq)) (denoted by asterixes 
connected by the full line) and the other vacuum expectation values introduced in sec- 
tion 5 as a function of the fraction / of molecules. The dashed line shows the square of 
the condensate. The expectation values of the 't Hooft operator ((0'*^°°^')) are repre- 
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sented by open squares. The behavior of the four quark condensates Oi — O4 is almost 
indistinguishable (full squares) and we show only one of them. All matrix elements are 
normalized to their value at / = 0. One clearly observes the restoration of chiral sym- 
metry as the fraction of molecules approaches / = 1. In our schematic model, not only 
chiral symmetry but also U{1)a symmetry is restored as / — 1. This is apparent from 
the behavior of ((0'*^°°^'^)), which is very large in the random model but vanishes in the 
molecular vacuumQ. The factorization assumption fails badly for this operator but works 
fairly well for the vector operators Oi. 

More microscopic information is given by the spectrum of the Dirac operator which we 
show in fig. 6. For the purely random vacuum, the spectrum is peaked at small eigenvalues 
A. Using the Casher Banks formula (qq) = —n{X = 0)/TT, this behavior shows that chiral 
symmetry is indeed broken. In the purely molecular phase, the distribution of eigenvalues 
is peaked at finite A, the density of eigenvalues at zero vanishes and chiral symmetry is 
restored. The position of the maximum of the eigenvalue distribution reflects the typical 
inverse size of the molecules. For arbitrary concentrations / we observe that the spectrum 
looks like a linear combination of a random and a molecular component. This means that 
the Dirac operator really decomposes into a random part and independent 2x2 blocks 
corresponding to molecules. 

Further details on the restoration of chiral symmetry are provided by hadronic corre- 
lation functions. Here we will only discuss spacelike correlators, since the corresponding 
screening masses have received a lot of attention in lattice calculations. We will present 
the corresponding temporal correlators in a forthcoming publication [Q. The correlators 
in the pion, delta, rho and ai channel are shown in fig. 7 and fig. 8. As the fraction of 
molecules approaches unity, one clearly observes that chiral symmetry is restored and the 
pion and delta as well as the rho and oi channels become degenerate. However, this does 
not mean that the correlation functions become purely perturbative. While the vector 
channels are indeed very close to the square of a massless thermal quark propagator, the 

^As explained in section 5, this does not preclude the possibility that the U{1)a symmetry remains 
broken in an unquenched calculation. 
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scalar meson correlation functions are significantly larger than what would be expected 
from the propagation of free quarks. 

This is also apparent from the results for the spacelike screening masses which we show 
in fig. 10. If the correlation functions are dominated by the propagation of free massless 
quarks, the spacelike screening masses are expected to behave as 27tT for mesons and SttT 
for baryons As the fraction of molecules approaches one, the screening masses for 
the vector mesons p and ai go to 27rT while the pion remains significantly lighter. This is 
consistent with the effective interaction discussed in the last section. Moreover, it is also 
in agreement with lattice results for the spacelike screening masses in the relevant regime 



T,~2T, 3 



We emphasize that the results shown in fig. 7 imply the restoration of U{1)a as well 



as chiral symmetry |39] since the isovector scalar (6) meson is not the chiral partner of the 



pion. As was remarked in this effect actually may have been observed on the lattice; 
the sigma meson correlation function determined in most lattice simulations is actually 
closer to the delta meson correlator. The reason is that the two loop graphs which 
distinguish the two correlation functions are difficult to measure and usually neglected in 
lattice simulations. 

In fig. 9 we show the behavior of the nucleon and delta correlation functions. Among 
the various possible correlation functions we only show the chiral even correlators of the 
first loffe current for the nucleon rji = e"''^(M"C7^M^)7^75(i'^ and the standard delta current 
jfi = e°''"^{u°'C'y^u^)u'^ (11^ and 11^ in the notation of |^). The chiral odd correlators 
vanish as / — 1, in agreement with chiral symmetry restoration. In this limit the chiral 
even correlators approach the cube of a thermal quark propagator. Similar to the pion 
channel, there is some residual interaction present in the nucleon channel. The screening 
mass of the delta approaches SvrT as / — 1, whereas the nucleon screening mass remains 
somewhat smaller. However, the effect is much less pronounced as compared to the pion 
channel. 

As a consequence of chiral symmetry restoration, the nucleon is expected to become 
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degenerate with its chiral partner, which is usually identified with the lightest odd parity 
N* resonance. Unfortunately, there is no current that couples exclusively to the A^* and 
not to the nucleon so that the corresponding screening mass is difficult to determine]^. 
However, the fact that the chiral odd nucleon correlation functions vanish implies that 
any contribution from a positive parity nucleon resonance to these correlators is exactly 
cancelled by the contribution of the corresponding negative parity resonance. This pro- 
vides at least an indirect check for the statement that the nucleon and its chiral partner 
are degenerate. 



7 Conclusions 

In this paper we have studied various consequences of a scenario in which the chiral phase 
transition in QCD is described as a rearrangement of the instanton liquid, going from a 
random phase (at low T) to a correlated phase of polarized // molecules at T > T^. In 
this scenario, a significant number of instantons is present at temperatures T = Tc ~ 2Tc, 
causing a variety of nonperturbative effects. 

As was shown in |13[ the phase transition does not require a perturbative instanton 



suppression but can be triggered by the temperature dependence of the fermion determi- 
nant. Here we have demonstrated that this particular behavior is due to the temperature 
dependence of the fermion overlap matrix elements. At temperatures T ~ these matrix 
elements strongly favor the appearance of polarized instanton— anti-instanton molecules. 
This correlated instanton configuration leads to chiral symmetry restoration. 

One important effect caused by instantons is the significant contribution to the energy 
density and pressure of the system near the phase transition region. As demonstrated in 
section 4 the rapid polarization of the instanton molecules causes the energy density to 
jump up, while the pressure lags behind and remains practically unchanged. This behavior 
is consistent with lattice simulations, in which the energy density rapidly reaches (or even 



^In lattice simulations one usually uses the fact that in staggered fermion calculations different parity 
components are discretized on alternating lattice sites in order to identify the odd parity component of 
the nucleon correlation function. 
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overshoots) the perturbative value whereas the pressure remains rather low until T ~ 2Tc 



25]. Thus, one can say that the energy of BNL AGS and CERN SPS accelerators is now 
used to a significant degree in order to produce gluomagnetic fields of non-perturbative 
origin! 

The presence of // molecules above Tc also produces quite specific interactions between 
light quarks. We have derived an effective local U{2)xU (2) symmetric fermion lagrangian 
which describes these effects. This Lagrangian is invariant under a remnant of the Pauli- 
Giirsey symmetry which results in the degeneracy of the rho and omega channels. We find 
that there is a substantial residual attraction in the scalar-pseudoscalar meson channels 
even above T^. This results is in agreement with lattice simulations, in which the presence 
of an attractive interaction in the scalar channel has been established from an analysis 



of spacelike screening masses |^ and direct measurements of the scalar susceptibility 
||4l[| . The considerably weaker interaction in the vector channel follows from the behavior 
of the quark number (vector) susceptibility |^2|, In agreement with our effective 



Lagrangian, lattice calculations result in a coupling constant that is four times larger in 



the scalar-pseudoscalar channel than in the axial-vector channel ^ 



We have also established the presence of molecule-induced interactions at T ~ Tc from 
direct measurements of the correlation functions in a schematic 'cocktail' model. In this 
model, the correlators are determined in an instanton ensemble with a given fraction / 
of correlated // pairs at temperatures around the critical temperature T ~ 150 MeV. As 
the instanton liquid becomes completely correlated, we observe the restoration of chiral 
symmetry. In this limit most of the spacelike screening masses approach their perturbative 
values, m = 2ttT for mesons and m = SvrT for baryons. Notable exceptions are the scalar 
(tt, a, 6, 7]') mesons whose screening masses are significantly smaller than 2iiT. 

Finally, we would like to comment that the 'mixed phase' (the region T ~ discussed 
above can be observed not only in lattice simulations, but also in ongoing real experiments. 
In fact, it is exactly the kind of matter which is now produced and studied at the CERN 
SPS and Brookhaven AGS (to be complemented by future experiments at RHIC and 

26 



LHC). The collective hydrodynamic expansion can naturally test the equation of state. 
Peaks in the dilepton spectra should provide information on the masses of vector mesons 
decaying in matter, so hopefully one can check modifications of the masses of the rho, the 
omega and the ai or the splitting of longitudinal and transverse components. 
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figure captions 



figure 1 Distribution of the polarization angle of an instanton— anti-instanton molecule 
for different temperatures T[MeV]. The distributions are shown as a function of the 
orientation angle a introduced in the text and are normalized to a uniform distribution 
cos^ OL. The two figures show the two cases of two light and one massive flavor (QCD) 
and four massless flavours {Nf — A). 

figure 2 Degree of polarization of an instanton anti-instanton molecule as a function of 
the temperature. The two figures show the two cases of two light and one massive flavor 
(QCD) and four massless flavors {Nf — 4). 

figure 3 Interaction and euclidean 'energy' < — > in units of the single instan- 
ton values for a molecule oriented in the timelike direction at a temperature T = 200 
MeV. Results are given as a function of the molecule size d in units of the instanton 
radius p = 0.35 fm. The dashed curves represents the results for the ratio ansatz whereas 
the full lines show the results for the sum ansatz. The lower line of each pair is for the 
most attractive orientation whereas the upper line is for the most repulsive orientation. 

figure 4 Interaction and euclidean 'energy' < — > in units of the single instan- 
ton values for a molecule oriented in the spacehke direction at finite temperature T = 200 
MeV. Results are given as a function of the molecule size d in units of the instanton radius 
p. Curves as in figure 3. 

figure 5 Quark condensates in the 'cocktail' model for various concentrations / = '^N^oi/Ninst 
with the total density Uinst — 0.75 fm~^ and the temperature T — 150 MeV fixed. The 
various operators shown are defined in section 5 of the text. All expectation values are 
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normalized to their value at / = 0. The solid and dashed lines shows the behavior of the 

quark condensate (gg) and the square of the condensate (qq)'^, the open squares represent 
the expectation value of the 't Hooft vertex, and the expectation values < Oi_4 > are 
denoted by the full squares. 

figure 6 Spectrum of the Dirac operator iD for various concentrations / =0.00, 0.50, 
0.75, 0.90, 1.00 of molecules. The density of instantons riinst — 0.75 fm~^ and the temper- 
ature T = 150 MeV are fixed. 

figure 7 Scalar meson correlation functions for various concentrations / =0.00, 0.50, 0.75, 
0.90, 1.00 of molecules. The density of instantons riinst — 0.75 fm~^ and the temperature 
T = 150 MeV are fixed. The dashed curve shows the correlation function for massless 
quarks. 

figure 8 Vector meson correlation functions for various concentrations / =0.00, 0.50, 0.75, 
0.90, 1.00 of molecules. The density of instantons riinst — 0.75 fm~^ and the temperature 
T = 150 MeV are fixed. The dashed curve shows the correlation function for massless 
quarks. 

figure 9 Nucleon and delta correlation functions for various concentrations / =0.00, 0.50, 
0.75, 0.90, 1.00 of molecules. The density of instantons riinst = 0.75 fm~^ and the tem- 
perature T — 150 MeV are fixed. The dashed curve shows the correlation function for 
massless quarks. 

figure 10 Spacelike screening masses as a function of the concentration / of molecules. 
The density of instantons riinst — 0.75 fm~^ and the temperature T — 150 MeV are fixed. 
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